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The effects of an automated, tenth order in k subtraction scheme on the noise variance of various Wilson QCD 
disconnected matrix elements are examined. It is found that there is a dramatic reduction in the variance of the 
lattice point-split electromagnetic currents and that this reduction persists at small quark mass. 



1. Introduction 

Many types of QCD matrix calculations require 
the global estimate of quark "disconnected" am- 
plitudes. Examples for the nucleon include the 
pi-N sigma term, quark spin content and the 
strangeness form factors. There is a continuing 
need for new methods which can extract these 
matrix elements more efficiently. 

Noise theory methods are based upon projec- 
tion of the signal using random noise vectors as 
input. That is, given 



Mx = rj, 



(1) 



where M is the N x N quark matrix, x is the 
solution vector and rj is the noise vector, with 



< rn >= 0, < rji^j >= 5ij, 



(2) 



where one is averaging over the noise vectors, any 
inverse matrix element, , can then be ob- 
tained from 



< VjXi >= M ik < W fe >= M y 

k 



— I 



(3) 



Subtraction methods can be of great assistance 
in reducing the noise variance Q. The key is using 
a perturbative expansion of the quark matrix as 
the subtraction matrices. The method described 
is completely iterative so that going to higher or- 
ders in perturbative subtraction is very easy and 
the overhead is also extremely low. 

2. The Method 

Let us review the results of noise theory for the 
expectation value and variance of matrices with 



various types of noises. Define § 
1 L 

Xmn = Y VmlVnl- 



1=1 

(to, n = 1, . . . , N; I = 1, . . . , L.) We have 

Y v* 

^■mn — ^ nm i 

and the expectation value, 



(4) 



(5) 



(6) 



The expectation value of Tr{QX} is Tr{Q} and 
it's noise variance is 

V[Tr{QX}} =< | ]T q mn X nm - Tr{Q}\ 2 > . (7) 



This results in 

V[Tr{QX}} =J2(< l X »™! 2 > l*> 



(8) 



nrn) + ^2 < \ X nn ~ 1| 2 > \ln 



(9) 



Consider three cases. First, general real noise: 

< |-X"mn| 2 > = < {X mn ) 2 >= — , 
for to ^ n so that 

V[Tr{QX rcal }} = jY, d«™»| 2 + ^q* nm ) (10) 

+ J2< - 1| 2 > knl 2 - 

n 

Compare this with the Z(2) case which also has 
Eq.(0) for to =/= n, but also 



< \X nn -l\ 2 >=0. 



(11) 
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This shows that 

V[Tr{QX m }] < V[Tr{QX real }] . (12) 
For the Z(N) (N > 3) case one has instead 
1 



< \X„ 



>= 



V 



< (x mn y >= o, 

for to ^ n, but again 

< \X nn - 1| 2 >= 0. 
Thus 



V[Tr{QX z(N) } 



| El* 



2 

mn | 5 



(13) 
(14) 

(15) 
(16) 



and the variance relationship of Z(2) and Z(N) 
is not fixed in general for arbitrary Q. However, 
if the phases of q mn and g* m are uncorrelated, 
then V[Tr{QX z{2) }} » F[Tr{QX zw }] (JV > 
3) , which is apparently the case for the operators 
studied here and in Ref. 
Now consider Q such that 



< Tr{Q} >= 0. 
Obviously, 

< Tr{(Q - Q)X} >=< Tr{Q} > 
However, 

V[Tr{(Q-Q)X}]^V[Tr{QX}}. 



(17) 



(18) 



(19) 



As we have seen for Z(N) (N > 2), the variance 
originates exclusively from off diagonal entries. 
So the trick is to try to find matrices Q which are 
traceless (or can be made so) but which mimick 
the off-diagonal part of Q as much as possible. 

The natural choice is simply to choose as Q the 
perturbative expansion of the quark matrix. This 
is given by ({I J} are collective indices) 



( M p 1 ){U} - j 5 1 

p 8{ij} - kP { ij} 

where 

x^y — a^ < 



(20) 



(21) 



Expanding this in k gives, 



M~ l = 5 + kP 



k 2 P 2 



k 3 P 3 



(22) 



One constructs < 7]j(M p 1 )ifc% > and sub- 
tracts it from < TjjM^Z T]k >, where r\ is the noise 
vector. This construction is an iterative process 
and so is easy to code and extend to higher powers 
on the computer. One can put coefficients in front 
of the various terms in Eq. ( p2[) and vary them to 
find the minimum in the variance, but such co- 
efficients take on their perturbative value || ex- 
cept perhaps for low order expansions^. Inter- 
estingly, significant variance improvement occurs 
in some operators even at 0th order in k. 

For a given operator, O, the matrix OM~ x 
encountered in < ^O^p > ga uge— Tr(OM~ 1 ) is 
not traceless in general. To correct for this one 
must re-add the perturbative part, subtracted 
earlier, to get the full, statistically unbiased an- 
swer. How does one calculate the perturbative 
part? Hard, exact way: explicitly construct all 
the gauge invariant paths (up to a given k order) 
for a given operator. Easy, statistical way: sub- 
ject the perturbative contribution to a separate 
Monte Carlo estimation, identical to the Monte 
Carlo applied to the nonpertubative part. This 
separate Monte Carlo is easy to do because one 
is simply constructing a matrix rather than in- 
verting one. Local operators require perturbative 
corrections starting at 4th order and point-split 
ones have corrections starting at 3rd order. In 
the following I will carry this procedure to k 10 . 

3. Variance Ratio Results 

I will show the ratio of unsubtracted noise 
variance to subtracted variance, V unsu b/V su b for 
Z(2) noise. Since computer time is proportional 
to the operator variance, the ratio gives a mea- 
sure of the decrease in the computer time needed 
to reach a given noise variance level. The lattices 
are Wilson 16 3 x 24, (3 = 6.0. Note that I am us- 
ing a "one noise" inversion method, meaning that 
a global estimate of the operator is obtained after 
a single matrix inversion. The appropriate oper- 
ators to apply this to are ipip and local and point 
split ipjuip (see Ref. ||). Fig. 1 shows the effect of 
the level of subtraction on the point-split charge 
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Figure 1. Ratio of unsubtracted to subtracted 
noise variance for zero-momentum point-split J4 
as a function of subtraction level at k = 0.148. 



Figure 2. Ratio of unsubtracted to subtracted 
variance for various operators at k — 0.148 af- 
ter 10th order subtraction. See the text for the 
meanings of the operators. 



density operator at n = 0.148. Fig. 2 shows the 
ratio of variances for the scalar ("S"), the four 
local vector ("L VEC 1-4"), and the four point- 
split vector operators ("P-S VEC 1-4") after 10th 
order subtraction, also at k = 0.148. 

4. Conclusions and Acknowledgments 

It has been demonstrated that a large reduction 
in the noise variance of certain lattice operators is 
obtained using perturbative subtraction methods. 
The method is effective for the scalar and local 
vector currents, but most effective for the point- 
split vector currents. The method can become 
less effective at lower quark masses, depending 
on the operator. The 10th order point-split vec- 
tor, local vector, and scalar variance ratios change 
from - 35, - 12, and ~ 10 at k = 0.148, to - 25 
<~ 10, and ~ 5 at n — 0.152, respectively. Similar 
methods can be devised for other operators (ax- 
ial, pseudoscalar, tensor) by implementing this 
algorithm in the context of "12 noise" methods. 

The operators condidered here are all zero mo- 
mentum. Of course for disconnected form fac- 
tor evaluations one is more interested in nonzero 



momentum data. Although the results are not 
shown here I have found essentially identical re- 
sults to the above for the momentum transformed 
data. These methods should be extremely useful 
in the lattice evaluations of strangeness form fac- 
tors, which are of current experimental interest. 

This work is supported in part by NSF Grant 
No. 9722073 and the National Center for Super- 
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gin 2000 System at the University of Illinois. The 
hospitality of the University of Kentucky Physics 
Department, where this work was begun, is grate- 
fully acknowledged. 
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